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Abstract 

In this note we present a study of the radiative tails in the invariant mass spectra 
of muon or electron pairs from J/ip,tp(2S) and T(IS') decays which is due to an 
additional emission of a photon. An analytic formula for dilepton mass spectra 
in radiative decays is derived and a Monte Carlo simulation for realistic detector 
conditions is used to study effects on the spectra. A rather simple parameterisation 
is given, suitable for the treatment of real data. 
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1 Introduction 



The prediction from the Standard Model for radiative decays of Z bosons into lepton 
pairs was discussed by Fleischer and Jegerlehner pQ and analytic results for the photon 
emission valid for arbitrary cuts on photon energies and angles have been presented. The 
lowest order decay width T (Z — > and contributions from radiative corrections were 
factorised in the presented formulae. We used the same formulae with appropriate mass 
parameters to estimate the radiative effects for heavy vector mesons decaying into lepton 
pairs. 

Radiative effects in the context of the bottomonia spectroscopy were discussed in 
details by Buchmiiller and Cooper [2j with a clear recommendation (p. 427) that the 
experimental definition of T — > e + e~ should include all T — > e + e~7 decays in order 
to avoid sensitivity to a cut in the photon energy. Certainly that should be valid also 
for charmonia decays and for fi + fi~ and / u + / u _ 7 final states, respectively. There are no 
principal differences between these cases, besides the level of the radiative effects. 

The first observation jH] of the radiative decay J/ip —>■ e + e~7, where the branching 
ratio B{J/tp — > e + e~j,E^ >100 MeV) divided by the branching ratio B{J/ip — > e + e~) 
was measured as 0.147±0.022, is in good agreement with the expectation. The theoretical 
prediction based on results pQ for the corresponding decay into muons is smaller by a 
factor about 3 only, i.e. also rather large. Radiative tails could be statistically significant 
in inclusive mass spectra of lepton pairs from decays of J/ ip in measurements with high 
statistics and rather good mass resolution like in the HERA-B experiment jlj. 

We derive an analytic formula for the dilepton invariant mass spectra in radiative 
decays. In practice, a convolution of such ideal spectrum with detector mass resolution 
should by taken into account. We investigated the smearing effects by Monte Carlo 
simulation taking into account realistic detector conditions for the J/ip, ip{2S) and T(1S) 
decays into /x + yU~7 or e + e~ r y. A simple parameterisation was found to describe simulated 
dilepton mass spectra for different mass resolutions. The parameterisation fitted rather 
well the simulated mass spectra and was appropriate for fitting of real dimuon spectra 
from J/ip decays. 

2 Photon Emission in Decays into Lepton Pairs 

In jT] radiative decays of Z bosons were considered and analytic results for arbitrary cuts 
on the photon energy and on the angles between photon and leptons were presented. The 
authors included the complete one-loop electroweak corrections and multi-soft-photon 
effects. The photons (real or virtual) could be emitted and absorbed only by the final state 
leptons. The effect of the initial state entered only via the phase space boundaries given 
by its mass, therefore only the mass of the decaying particle entered into the expressions 
for radiative corrections and the lowest order decay width T (Z — > was used as 

a normalisation factor. We used the appropriate formulas from pQ to study radiative 
decays of heavy vector mesons into lepton pairs. 
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Let us consider the decay X — > l + l of a heavy vector state with mass M into a pair 
of leptons with masses m/. The lowest order decay width is 

r = r (x -> /+/-). (i) 

The bremsstrahlung process 

^(Po)-r(Pi)+/ + (p 2 )+7W (2) 
is distributed in the phase space as PQ: 
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where a is the fine structure constant, c = {j>\ + P2) 2 /M 2 , r = (p — pi) 2 /M 2 and 
a = Amf/M 2 , a < q < 1. 

In the rest frame of the decaying particle we have q = 1 — 2E 7 /M, r = 1 — 2Ei/M 
and 1 — ^ — r = 1 — 2E 2 /M, and in the center of mass system of the lepton pair: 

r = 1 — s— -(1 — <?)(! — rcosfli) or r = ^(1 — <?)(! - r cos0 2 ) 



with r = y 1 — a/q and 6*j the angle between the photon and lepton with momentum p i: 
obviously 6\ = n — 6 2 - Variations of 62 from 62 = to 6 2 = it corresponds to 

r(0) = ^(l-0(l-r), r(n)= 1 -(l-q)(l + r). 

The parameter <j is related to the invariant mass of the lepton pair: 

m 2 = ( Pl +p 2 ) 2 = q- M 2 . 

The distribution P(q) can be evaluated by integrating the distribution P(q,r) over the 
parameter r: 

1 dT(X -> Z+Z-7) , /-W _ , , 



T dq J T (o) 

a 1 + q 2 /, 1 + r 



, In r , (4) 

The latter formula agrees with the expression P(q, > 5) from jTj in the limit 5 = 
presenting the result of integration of P(q,r) over the angles 61,62 with exclusion of 
cones around the leptons defined by the angle 5. 

Photons with sufficient energy can be detected. The fraction of decays corresponding 
to the emission of hard photons is 

1 r l-2E min /M 

C hard (E min ) = — ■ T(X -> £ 7 > E mm ) = / Pis) dq, (5) 

1 J a 



where E min is the minimal photon energy in the rest frame of the decaying particle. The 
result of integration for E min «C M/2, neglecting terms of order 0(aa), is known [I]: 
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The measurement of — > e + e 7 [3] is in good agreement with the calculation 
in fig. [TJ 
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Figure 1: Parameter Chard as a function of the minimal energy E m i n of the photon in the J/ip 
rest frame for radiative decays J /if) — > /x + /i~7 (solid line) and J/ip — > e + e~7 (dashed line). The 
point with errors was evaluated from the E760 result . 



The lepton mass entered in eq. |H1 in the first (leading) term as ln(M 2 /m 2 ) which is 
smaller only by a factor 2.6 for muons compared to electrons. The effect of hard photon 
emission is expected to be rather large also for decays into yU + /i~ 7 as shown in fig. ^ The 
effects are larger for heavier decaying particles as seen in fig. El comparing J/ip, ip{2S) 
and T(IS'), respectively. 

In order to get the physical rate, we have to take into account contributions from 
virtual photons (emitted and absorbed by leptons) and emitted soft photons (E 1 < E min ). 
The sum of these contributions is known pQ: 

C * soft\Erniri) 

The total width taking into account / + /~, 7 configurations is 

r M (x -> i + r, z + r 7 ) = (i + c soft + c hard ) ■ r (x -> i+r). (8) 

The term C so f t is negative as seen in fig. El The leading terms in C so f t and Chard have 
singular behavior for E min 0, but in the sum C so f t + Chard they exactly cancel, as it 
should be: 

Cgoft + Chard = ^ 77- (9) 
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Figure 2: The parameter Chard as a function of the minimal energy E m i n of the photon in 
the rest frame of J ftp (solid line), ip(2S) (dashed line), and T(15) (dotted line) decaying into 
/i + ^ _ 7 (a) or e + e~7 (b). 
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Figure 3: The parameters Chard (solid line) and C so ft (dashed line) for the radiative decay 
J/i/j — > / u + / u _ 7 as functions of the minimal photon energy E m i n in the rest frame of the J /ip. 



When an experiment measures all / + /~, 7 configurations independent of the en- 
ergy of emitted photons, the effects of radiative corrections appear as a small increase of 
the decay width 

= -- « 0.00174. (10) 
To 2n 2 

However, cuts on the energy of the emitted photon can have much stronger effect. Let us 
assume that we are selecting lepton pairs inclusively and apply a cut on their invariant 
mass 

\m-M\ < A. 
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The radiative decays with photon energies in the rest frame of the decay 



^ >Ba = A ( 1 -4) 

will be rejected and the number of selected decays will be reduced by the factor 

N(l+r, \m-M\ < A) (l + C soft (E A )) 



'IV. 



N(l+l~) 



(l + f 1 



(12) 



For inclusive registration of fi + fi~ pairs from J/ip decays and the rather soft cut A = 
200 MeV, about 3% of the decays will be rejected. In general, the dilepton mass m is 
shifted by photon emission 



m 



>M(M - 2£ 7 ) « M - £ 7 (£ 7 <M). (13) 
Hard photon emission causes a tail towards lower masses in the dilepton mass spectrum. 



10' 



io- 



10 



10 



10 




0.5 1 1.5 2 2.5 3 
Dimuon mass [GeV] 



Figure 4: The dimuon mass spectrum defined by eq. El (solid line) and results of Monte Carlo 
simulation (points) of the decay J/ip —* for an ideal detector (without smearing, 4-7T 

acceptance) for E m i n = 10 MeV. 



The distribution P{m) of the dilepton mass in the radiative decay (J2J we find from eq. (jU 

as 



i dv{x ^m~i) 

r n dm 



P(m] 



a 2m 



m A \ / 1 + r 
7r (M 2 — m 2 ) M 1 / V n I^7~ r 



(14) 



where r = y 1 — 4m; 2 /m 2 is also a function of m. The singular behavior of the spectrum 
for E 7 — > (m — > M) is clearly seen in fig. 0] 
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In an experiment which measures the dilepton mass spectrum inclusively a separation 
of and / + /~7 final states is not possible on an event by event basis. For very low 
photon energies that is not possible even on a statistical basis, because both decays 
contribute to the peak around the mass of the decaying particle. The measured spectrum 
can be described by two components. The main component corresponds to decays into 
l + E and / + Z~7 (E 1 < E min ) contributing to a "monochromatic" line (m = M) in dilepton 
mass spectrum smeared by detector resolutions. An additional component is related to 
the / + / _ 7 (2£ 7 > E min ) final states producing the spectrum dHj) also affected by detector 
resolutions. The fraction of events in the latter component is 

W hard = \ 3 • (15) 

' 2tt 2 

3 Smearing Effects due to Detector Resolutions 

After smearing due to limited detector resolutions, the shape of the dilepton mass spec- 
trum from radiative decay (|14|) changes in particular for low energy photons (m rs M) as 
seen in fig. The spectrum (fT4"j) convoluted with the Gaussian could not be expressed in 
an analytic form and in addition, experimental cuts complicate the situation even more. 
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Figure 5: Dimuon mass spectra in the radiative decay J ftp — ► 7 for E m i n = 10 MeV 

(solid line) smeared with mass resolution of about 42 MeV (points connected by dashed line). 



Our goal was to derive a function to be used in fitting the spectra. We investigated 
radiative decays by Monte Carlo simulation for real experimental conditions and propose 
a rather simple parameterisation for smeared dilepton mass spectra. We used the Monte 
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Carlo generator for simulation of the kinematics of the decay (J2J) developed by Lohse [5] 
and based on the results in pQ. The generated spectrum of the dilepton invariant mass 
was in good agreement with the analytic expression (jl4j) as shown in fig. |3J To study 
the real experimental conditions we simulated radiative decay within the framework of 
the HERA-B detector. Protons with an energy of 920 GeV colliding with a nuclear 
target directly produced the J/ip (or ip(2S), T(15)). The distributions of the transverse 
momentum p? and the Feynman scaling variable xp were generated according to the 
model described in [5]. The particle decayed into n + n~ (e + e~) or fi~ r y(e + e~ r y) with a 
probability determined by the cut parameter E m i n , the minimal energy of the photon in 
the rest frame of the decay. 

The HERA-B detector is a forward spectrometer with a tracking system including 
the vertex detector placed in front of the magnet and main tracker system behind the 
magnet. The electro- magnetic calorimeter and muon system located behind the main 
tracker are used for particle identification and triggering. The main tracker includes the 
inner and outer trackers, but only the latter was taken into account since only this part 
was used in the dilepton trigger. Events where the both leptons passed through the 
tracking and particle identification detectors were considered. The trigger also imposed 
cuts on leptons, e.g. only muons with energy higher then 6 GeV were selected. 

The momentum resolution for muons was simulated according to the parameterisa- 
tion [7j 

— (%) = 0.94 © 9.55 • 10~ 3 • p © 0.484 • p a263 , (16) 
p 

where p is in GeV. We were interested to isolate effects induced by the internal bremsstrahlung 
in decays and therefore did not simulate bremsstrahlung photons from electrons on the 
path through the detector material. The simulated momentum resolution for electrons 
was the same as for muons. To study dilepton spectra for various resolutions we multi- 
plied eq. Uniby fudge factors from 0.7 up to 2.5. That corresponded to the variation of 
mass resolution at the J/ip peak from 23 MeV up to 81 MeV, respectively. 

The decays corresponding to E 1 < E min were simulated as final states smeared 
by detector resolutions. The remaining part of the spectrum consisted of simulated Z + /~7 
configurations (E 1 > E min ), also affected by detector resolutions. Our aim was to find 
a simple analytic approximation for the latter part of the simulated spectrum with the 
additional condition that the mass resolution can also be treated as unknown and included 
in the fit parameters. We solved this problem in three steps: 

• Find the value E min appropriate for the simulation; 

• Define a parameterisation for the smeared dilepton mass spectra in the radiative 

decays; 

• Determine the corresponding parameters as functions of the mass resolution. 
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Figure 6: Dimuon mass spectra simulated by Monte Carlo for the mixture of J/tp — ► 
and /i + /i~7 decays. The latter decay was evaluated with the probabilities corresponding to 
Emin =10 MeV (points) and 200 MeV (histogram), respectively. The simulated mass resolution 
was about 42 MeV at the J/ip peak. 
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Figure 7: Fraction of events in the tail of dimuon mass spectra lower then 100 MeV as the 
J/tj} mass as a function of the cut parameter E m i n . This fraction was evaluated for the mixture 
of J/i/j — > fi + fi~ and / u + /x _ 7 decays simulated with probabilities corresponding to E m i n . The 
mass resolution was about 42 MeV at the J ftp peak. 
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On the one hand, E min should not be extraordinarily small to stay away from the 
singularities in eqs. ()6l7|) . On the other hand, the simulation with E min much bigger than 
the mass resolution produces unreasonable mass spectra, as is clearly seen in fig. El for 
the spectrum simulated with E m i n = 200 MeV and a mass resolution of 42 MeV. If the 
photon energy is small enough compared to the mass resolution a, it does not matter 
whether or Z + Z~7 was simulated, because the invariant mass is smeared with the 
resolution a ^> M — m. The shape of the inclusive dilepton mass spectrum becomes 
stable against a decrease of E min . 

To make a proper choice we estimated the fraction of events in the "tail" of the dilep- 
ton mass spectra (m < M(J/ip) — 100 MeV) as a function of the parameter E min for 
decays J/ij) — > /i + /i~7 simulated with the mass resolution 42 MeV. This depen- 

dence, presented in fig. shows a plateau where the results are stable with respect to 
the cut parameter E min . We selected the value E min = 10 MeV for further simulation. 



4 Parameterisation of Smeared Mass Spectra 

In this section the parameterisation of the smeared dilepton mass spectrum from the 
radiative decay (J2J) with E^ > E min is given. The upper edge of the spectrum is the 
result of smearing of the sharp peak in the dilepton mass spectrum ()14j) atmfsM and 
it looks like a Gaussian (see fig. EJ). For the full dilepton mass range we propose the 
following parameterisation: 

V{rn) = G(m\mR, or), m > ms (17) 

and 

V(m) = cs ■ G(m\mR, asijn)), m<ms, (18) 
where G(m\mR, ctr) is the Gaussian function 

nf | \ 1 ( {m-m R ) 2 \ 

G{m\m R , a R ) = —= — exp — 

V2ticxr \ 2a r J 

with an average tur and standard deviation aR, and G(m\mR, as{m)) that with standard 
deviation 05 (m), respectively, and cs is a scaling factor. For a given mass resolution, the 
parameters rriR,aR and C5 are considered as constants, but 0"s(m) as a polynomial: 

a s (m) = Y j a 2 i {m R -m) i (19) 

i 

with constants a« to be defined. Each term in the polynomial is positive which helps to 
reduce mutual correlations of parameters during their search. 

In the point m = ms we require continuity of the function V(m) and its first derivative 
Vim). From the first condition we define the parameter C5: 

&s(ms) ({rn s -m R ) 2 (m s - m R ) 2 \ 

°s = exp 2( — : — -2 • (20) 

a R \ 2a 2 s {m s ) 2a 2 R J 
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From the second condition follows the nonlinear equation 

1 das{ m ) m — m R (m — mn) 2 das(m) 
crs(m) dm a s( m ) a s( m ) dm 

to define mg. 
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Figure 8: Monte Carlo results for the dimuon mass spectrum from the radiative decay J /ip — > 
/U + /*t~7 (histogram) and the parameterisation described in text (solid line). The simulation was 
made for a mass resolution of about 42 MeV at the J ftp peak. 

The remaining parameters mjj, <7r, Oj were used as free parameters and determined 
from the condition of the best fit of V(m) to the dilepton mass distribution for 10 6 
simulated radiative decays. For the decay J /if) — > the parameters 01,03,05,010 

were found to be statistically significant. The function Vim) with these parameters 
describe the simulated dimuon spectrum quite well, as seen in fig. |H1 

We repeat the fitting of spectra simulated with different detector momentum reso- 
lutions and evaluated the dependences of the considered parameters as functions of the 
mass resolution a = a(m) for m = M. Linear or quadratic polynomials in a were suffi- 
cient to define these functions. The parameter a is determined from the Gaussian fit of 
the /x 4 'f/T mass spectra for independently simulated J/ip — > fi + fi~ decays. 

The parameter <jr depends linearly on the mass resolution as shown in fig. UJi). The 
nonlinearity of the parameter mjj is quite small as seen in fig. ^jjp). The point 777.5, 
where the Gaussian distribution G(m, thr, or) turns to the long tail, moves towards 
lower masses for a worse mass resolution (see fig.HOb). The coefficients of the polynomial 
([19)1 were rather well approximated by quadratic functions as presented in fig. HUb) for 
the parameter oi. 
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Figure 9: The parameters <jr (a) and rriR (b) for the dimuon spectrum from the decay 
J/ip — > / u + / u _ 7 as functions of the mass resolution <r. 
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Figure 10: The parameters mj (a) and a\ (b) for the dimuon spectrum from the decay 
J/ijj — > / u + / u~7 as functions of the mass resolution <r. 



Finally, we approximated as a quadratic polynomial in a the factor 

1 



N R (a) = 



J V(m\a) dm 

to be used for normalisation of the function V(m\a) for the given mass resolution. 



(22) 
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As an example, we present the functional dependence for parameters determined for 
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where a is varying from 0.023 to 0.081 GeV. 
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Figure 11: Monte Carlo results for the dimuon mass spectrum from the radiative decay 
T(IS') — > (histogram) and their parameterisation (solid line). The simulation was 

made for the mass resolution of about 176 MeV at the T peak. 



The parameterisation describes also well the dimuon mass spectra simulated for the 
radiative decay T — > // + /^~7 as seen in fig. ^2 The shape of the tail looks more com- 
plicated than for the J/ip. The shoulder around 4 GeV can be explained by acceptance 
effects for the simulated detector geometry. These effects appear in the mass range cor- 
responding to very low values of V(m\a). An additional parameter a 2 o was needed in 
polynomial ()19|) to overcome the complications of the mass spectrum. 



13 



Table 1: Sets of coefficients of the polynomial (fHI)) found to be sufficient to parameterise 
inclusive dilepton mass spectra simulated for radiative decays of different particles. 



Particle 


Decay into fi + fi 7 


Decay into e + e 7 




01, 03) fl 5) 010 


0i ; 03, 05, aio, ct2o 


ij(2S) 


01) d3; fl 5) 010 


0i ; 03, 05, aio, ct20 


T{1S) 


01) 03) 05) 010) 020 


01) GS3, a 5) a 10, «28 



We simulated radiative decays of the J/i(),if)(2S), T(1S) into muon or electron pairs 
and evaluated the corresponding parameterisations for the dilepton spectra. In all cases, 
the parameters a%, 03, 05, a w were found to be significant as shown in table [T] For decays 
with smaller rrii/M an additional parameter was needed (020 or g^s), to take into account 
specific details induced by the simulated experimental set-up. These parameters are 
needed at the lower edge of the mass spectra where the probability densities are low. 

5 Treatment of Inclusive Measurements of Dilepton 
Mass Spectra 

The dilepton mass spectra measured inclusively can be described as 

5(771 \M,a) = (l- W hard ) ■ G(m\M, a) + W hard ■ N R (a) ■ V(m\ M, a), (23) 

where V(m \ M, a) is the function (| 17118)) with corresponding parameterisation normalised 
by the factor N R (a) and the weight Whard defines the relative contribution of the radiative 
tail into the spectrum. The first term describes the configurations without photon and 
with low energies photons emissions smeared accordingly by detector resolutions. In a 
real experiment this term could be a more complicated function than a Gaussian. For 
electrons also an additional term could be appropriate to describe the bremsstrahlung 
process on the path through the detector material. As it was mentioned already, this 
process was not simulated for the present study. The weight Whard can be defined as 

Chardi^Emin) ' ■^■hard\-^min) (24) 

(1 ('sof t ( ^Vrmi) ) ' J ^soft\-E^min) Chard^J^min) ' ^hardi^min) 

where A so f t (E min ) is an acceptance for l + l~ and (E 7 < E min ) final states and 

A h ard(Emin) that for / + Z~7, > E min ). For close A soft and A hard the definition of 
Whard by eq. (fT5j) can be sufficient. The latter definition for E min = 10 MeV we use for 
further calculations. 

The corresponding mass spectra for the J/ip decays are presented in fig. El The inte- 
gral of the function ()23|) over the full interval of dilepton mass equals unity by definition. 
For the fitting of the histogram we used the product of the corresponding function 
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Figure 12: Monte Carlo results for inclusive dilepton mass spectra from decays of J/ip into 
e + e _ , e + e~7 (histogram) and /i + yT , /j> + fi~j (points) and the fit curves using the parameter- 
isations (|23|) fsolid curves). The mass resolution was about 42 MeV at the J/tp peak for both 
decays. 

and the number of events treated as fit parameter. The function (f2*3*J) itself depends on 
two parameters: the mass average M and resolution a of the Gaussian part. Hence, there 
are three fit parameters in total. 

The fitted curves describe well the results of the simulation as shown in fig. ED The 
fitting procedures for decays of the ip(2S) shown in fig. IT3k) and T(IS') in fig. CHb) were 
analogous to those described above. The fitting function also agreed with the simulation 
in a wide range of invariant masses. 

To make a test on real data, we used the HERA-B data recorded with a dimuon 
trigger during the 2002-2003 HERA running period. Figure El shows the dimuon mass 
spectrum obtained from these data with a prominent signal above some background. The 
background was described by an exponential function. 

In the HERA-B detector, a forward spectrometer with non-negligible amount of ma- 
terial on the path through the magnet jH], Moliere scattering is essential f?] for the 
momentum resolution. Therefore in the fitting function (}2*3*|) instead of one Gaussian we 
used a symmetric function composed of three Gaussians to take into account the details 
of the momentum (mass) resolution in the HERA-B detector. The fit of the dimuon 
mass spectrum by this symmetric function accompanied by an exponential function for 
the background resulted in \ 2 — 338 for 108 degrees of freedom as presented in figIT4*k). 
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Figure 13: Monte Carlo results for inclusive dilepton mass spectra from onia decays into 
e + e~, e + e~7 (histogram) and /x + /i~7 (points) and parameterisations (|23j) (solid curves). 

a) Decays of the ip(2S) were simulated with a mass resolution of 52 MeV at the signal peak. 

b) Decays of the T(1S) for a mass resolution of 176 MeV at the signal peak. 
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Figure 14: The invariant mass distribution of /x + around the J/ip peak measured in the 
HERA-B experiment and fitted without the radiative tail taken into account (a) and with 
the parameterisation including the radiative tail (b). The values of x 2 were 338 and 112, 
respectively, for 108 degrees of freedom. The background was described as exponential function 
(dashed curves). In inserts are shown: number of entries in the histogram, x 2 , fitted number of 
J ftp, fitted mass and mass resolution in GeV. 
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The inclusion of the radiative tail into the fitting function improved the result as seen 
in figJ14b). where \ 2 — 112 for 108 degrees of freedom. One should remember that the 
radiative tail (second term in (|23jl) was defined through parameters of the first term, i.e. 
M and a of the main peak, and the number of fitted parameters, and thus the degrees of 
freedom, were the same. Taking the radiative tail into account in the inclusive spectra of 
the dilepton mass could be significant for large statistics and rather high mass resolution. 

6 Summary 

Well known analytic results for the radiative decay of the vector boson Z — > / + Z _ 7 valid 
for arbitrary cuts in the photon energy and angles were used to estimate contributions of 
similar decays for charmonia and bottomonia vector states. The contributions were found 
to be rather large not only for decays into electrons but also into muons. The radiative 
effects appeared as tails in inclusive mass spectra of the lepton pairs. An analytic for- 
mula for the dilepton mass spectra in radiative decays is derived. The spectra convoluted 
with the detector mass resolution could be described by rather simple parameterisations. 
The parameterisations were studied by Monte Carlo for J/tp,ip(2S),T(lS) radiatively 
decaying into /i + /i~7 or e + e~7 using a simulation describing the real experiment. The 
parameterisation described the Monte Carlo results quite well and was found to be ap- 
propriate also for fitting real data. 
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